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Abstract 

A parton shower method for the photonic radiative correction is 
apphed to the single VF-boson production processes. The energy scale 
for the evolution of the parton shower is determined so that the correct 
soft-photon emission is reproduced. Photon spectra radiated from the 
partons are compared with those from the exact matrix elements, and 
show a good agreement. Possible errors due to a inappropriate energy- 
scale selection or due to the ambiguity of energy scale determination 
are also discussed, particularly for the measurements on triple gauge- 
couplings. 
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1 Introduction 



Single-iy production processes present an opportunity to study the anoma- 
lous triple gauge-couplings (TGC)llI] in the experiments at LEP2 and at fu- 
ture e^e~ linear colliders. For their precise predictions of cross sections 
the inclusion of an initial state radiative correction(ISR) is inevitable in the 
event-generator. As a tool for the ISR the structure function(SF)0 and 
the parton shower^ methods has been widely used for e^e~ annihilation 
processes. For the case of the single-VF production processes, however, the 
main contribution comes from non-annihilation type diagrams. The uni- 
versal factorization method used for the annihilation processes is obviously 
inappropriate. The main problem lies in choosing the energy scale of the 
factorization. A previous study on the two photon process @] has shown that 
the SF and QED parton shower(QEDPS) methods are able to reproduce pre- 
cisely the exact 0{a) results even for non-annihilation processes, as long as 
the appropriate energy scale is used. 

In this report a general method to find the energy scale for SF and QEDPS 
is proposed. Then numerical results of testing SF and QEDPS for e~e+ 
e'Veud and e'e"*" — > e~z>e/i'^z/^ are presented. Systematic errors are also 
discussed. 

2 Calculation method 

2.1 Energy Scale Determination 

The factorization theorem for QED radiative corrections in the leading- 
logarithmic approximation is valid independent of the structure of the matrix 
element of the kernel process. Hence SF and QEDPS must be applicable to 
any e^e~ scattering processes. However, choice of the energy scale in SF and 
QEDPS is not a trivial matter |1TT|. For simple processes considered so far 
like e^e~ annihilation and two-photon process with only multi-peripheral di- 
agrams, the evolution energy scale could be found by the exact perturbative 
calculations. However, this is not always the case when more complicated 
processes are concerned. Hence a way to find a suitable energy scale without 
knowing exact loop calculations should be established somehow. 

First let us look at the general consequence of the soft photon approxi- 
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mation. The cross section with radiations in the soft-photon hmit is given by 
the Born cross section multiphed by the following factor up to the double-log 
term as 
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where mj{pj) are the mass (momentum) of the j-th charged particle, kc the 
maximum energy of the soft photon (the value to separate soft- and hard- 
photons), E the beam energy and ej the electric charge in unit of the e"*" 
charge. The factor r]j is —1 for the initial particles and -|-1 for the final 
particles. The indices run over all the charged particles in the initial 
and final states. 

For the two-photon process, e~(]?_)e"'"(p+) e~(q'_)e+(g+)/i~(A;_)/i+(A;+), 
it is shown in ref.Q] that the soft-photon factor in Eq.(0) with a (p_ • g_)- 
term reproduces the 0(a) corrections^ up to the double-log term in the 
soft-photon limit. This implies that one is able to read off the possible evolu- 
tion energy scale in SF from Eq.([l|) without doing explicit loop calculations.^] 
The point is the observation that the energy scale s = (p_ + p+)^ does not 
appear in the soft-photon corrections even they are included in the general 
formula Eq.(|ip. In the case of the two-photon process we have ignored those 
terms in SF which come from the photon bridged between different charged 
lines, because the contributions from the box diagrams with a photon ex- 
change between e"*" and e~ is known to be small0. Fortunately the infrared 
part of the loop corrections is already included in Eq.(||) and there is no need 
to know the full form of the loop diagram. Let us look at two terms with, for 
example {p- ■ p+)- and (g_ • ]?+)-terms. The momentum of e~ is almost the 
same before and after the scattering(p_ ~ g-). Only the difference appears 
in Tjjrjk = +1 for a (p_p+)-term and r]jr]k = —1 for a (g_p+)-term. Then 
these terms compensate each other after summing them up for the forward 
scattering which is the dominant kinematical region of this process. This is 

^ A similar idea is independently proposed by Montagna et al. inlsl. 
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the mechanism that the energy scale s = (p_ does not appear in the 

soft-photon correction despite that it exists in Eq.(|l]). 

When some experimental cuts are imposed, for example the final e~ is 
tagged in a large angle, this cancellation is not perfect but partial and the 
energy scale s must appear in the soft-photon correction. In this case the 
annihilation type diagrams will also give a contribution. Then it may happen 
that the usual SF and QEDPS for the annihilation processes are justified to 
be used for the ISR with the energy scale s. To find the most dominant 
energy scale under the given experimental cuts an easy way is to integrate 
numerically the soft-photon cross section given by Eq.(|I]) over the allowed 
kinematical region. Thus in order to determine the energy scale it is sufficient 
to know the infrared behavior of the radiative process using the soft-photon 
factor. 

Let us determine the energy-scale of the QED radiative corrections to the 
single-iy production process, 

+ e+(p+) e'{q^) + i'e{qu) + u{ku)+d{kd). (3) 

The soft-photon correction factor in Eq.([^) is numerically integrated with 
the Born matrix element of the process only with the i-channel diagrams 
without any cut on the final fermions. In order to separate the contribution 
from each term, we take terms up to 0{a) in Taylor expansion of exponen- 
tial function in Eq.([^). Parameters used in the calculation will be explained 
in section 4. The results are shown in Table |l[ One can see that the main 
contribution comes from an electron charged-line (p_g_-term) and a positron 
charged-line (p+/CuA;rf-terms=p+A;„-term -|- p^kd-teim + fc^fc^-term), while all 
other contributions are negligiblly small. Like the two-photon processes the 
energy scale s does not appear in the soft-photon correction. The above 



all terms 


P-(l- 


P+kukd 


all other combinations 


1 


0.38 


0.61 


1.9 X 10"^ 



Table 1: Soft-photon correction factor from sets of the charged particle combinations for 
the process of e+e" — > e^Veud at the CM energy of 200 GeV. The factor from all terms is 
normalized to mrity. kc = IGeV is used, p+fcufc^^p+fcu-term + p+fcrf-term + fcufc^-term. 

results clearly indicate that one should apply SF or QEDPS to the elec- 
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tron and positron charged- lines individually with the energy scale of their 
momentum-transfer squared. 



2.2 Structure Function Method 

The analytic solution of the DGLAP evolution function [|1^] in the leading- 
logarithmic order is known as the structure function]^. With SF the QED 
corrected cross section is given by 

cTtotaiis) = J dxi^ J dxp^ J dxj+ J dxu J dxdDe~{xj.,-t^)De~{xF-,-t-) 
De+{xi+, ■~t+)Du{xu, Sud)Ddixd, Sud)cro{s), (4) 

where ctq is the Born cross section and D's are the SF. The energy scales 
t- = {p- - t+ = {p+ - {ku + kd)f and s^d = {K + kdf are chosen 
following the result of the last section. After(before) the photon radiation 
the initial(final) momenta become p± {q±) 

p- = xi-p-, q_ = ^—q_,... (5) 

Xp^ 

respectively. Then the CM energy squared, s, is scaled cLS S — Xj—Xj-\-S. 

2.3 Parton Shower Method 

Instead of the analytic formula of SF a Monte Cairo method based on the 
parton shower algorithm in QED can be used to solve DGLAP equation in the 
LL approximation. Its detailed algorithm of the QEDPS is found in ref.|T2 



for the e^e~ annihilation processes, in ref.|]T3| for the Bhabha process, and in 
ref.0 for the two-photon process. The same energy scale as the SF method 
is used in QEDPS also. Significant difference between SF and QEDPS is 
that the QEDPS can treat the transverse momentum of the emitted photons 
correctly by imposing the exact kinematics at the e — > 67 splitting. It does 
not affect the total cross sections so much when the final e~ has no cut. 
However, the finite recoiling of the final can result some effects on the 
tagged cross sections. 

In return for the exact kinematics at the e ^ 67 splitting, are no more 
on-shell after a photon emission. On the other hand the matrix element of the 
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hard scattering process must be calculated with on-shell external particles. 
A trick to map the off-shell four-momenta of the initial to those at on-shell 
is needed. The following method is used in the calculations. 

1. s = {p- + pj^y is calculated, where p± are the four- momenta of the 
initial after the photon emission by QEDPS. s is positive even for 
the off-shell e^. 

2. New four-momenta of initial in their rest-frame, p±, are defined as 
p^ = ml (on-shell) and s = (]?_ + All four- momenta of final 
particles are generated in the rest-frame of p+ + p_ . 

3. All four-momenta are rotated and boosted to match the three-momenta 
of p± with those of p±. 

This method respects the direction of the final rather than the CM energy 
of the collision. The total energy does not conserve then because of the 
virtuarity of the initial e^. The violation of the energy conservation is of the 
order of 10^^ GeV or less, and the probability of the violation more than 1 
MeV is 10-^ 



3 Numerical Calculations 



3.1 Cross sections with no cut 



Total and differential cross sections of semi-leptonic process e~e'^ e~i>eud 
and leptonic one e~e+ e~h'efJ-'^i^fi, are calculated with the radiative cor- 
rection by SF or QEDPS. Feynman diagrams of the semi-leptonic process 
are shown in Fig.^. Fortran codes to calculate the amplitudes are automat- 



ically produced by GRACE system[14|. All fermion-masses are kept finite in 



the calculations. Numerical integrations of the matrix element squared in 



the four-body phase space are done using BASES ||I5[. For a test with no ex- 
perimental cuts, only the t-channel diagrams (non- annihilation diagrams) are 
taken into account. Standard model parameters used in the calculations are 
summarized in Table^. The on-shell relation strict at the tree-level is em- 
ployed to determine weak couplings. Some of electro-weak corrections could 
be taken into account through the G^-scheme and running coupling constant. 
In this report, however, we do not include those effects because here we are 
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interested in looking at the pure QED radiative corrections to the Born cross 
section. It has been pointed out that a tiny violation of the gauge invariance 
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80.35 GeV 


Tw 


Mz 


91.1867 GeV 
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a 


1/137.0359895 


sin^e 


me 


0.511 X 10-3 GeV 




m„ 


5.0 X 10-3 GeV 


rrid 



1.96708 GeV 
2.49471 GeV 

105.658389 x 10"^ GeV 
10.0 X 10-3 GeV 



Table 2: Standard-model parameters 

caused by the inclusion of a finite VF-boson width results a wrong cross sec- 
tion for single- production |l^, |T^. To cure this problem the fermion-loop 



scheme I T7, 18, 19 1 has been proposed. It is reported in ref.[18|, however, 
that no significant difference is seen between the fermion-loop scheme and 
the fixed width scheme. Thus in this report the fixed width scheme, together 



with an appropriately modified current, described in |jT6|] is employed. 

For the total energy of the emitted photons both methods, SF and QEDPS, 
must give the same spectrum once the same energy scale is used. This is con- 
firmed for the semi-leptonic process as shown in Fig.H at the CM energy of 
200 GeV. We choose the same energy-evolution scale as SF described by 
Eq.(^). Total cross sections as a function of the CM energies at LEP2 with 
no cuts are shown in Fig.^. The effect of the QED radiative corrections on 
the total cross sections are obtained to be 7 to 8% in the LEP2 energy re- 
gion. If one used an inappropriate energy scale, say s, in SF, the ISR effect 
is overestimated about 4% as seen in Fig.|^. The result of SF is consistent 
with the QEDPS around 0.2% if the proper energy scale is employed. 

There is some ambiguity to choose the energy-evolution scale in the 
leading-log order. For example the transverse-momentum squared or the 
invariant-mass squared of ud (or /iz/^) systems could be a candidate of the 
energy scale. It is found that even these energy scales are chosen in QEDPS 
instead of {p+ — {ku + kd)Y, the total cross sections change by only 0.6%, 
which must be allocated to the theoretical uncertainty. 

The energy and angular distributions of the hard photon from QEDPS 
are compared with those obtained from the exact matrix elements. The 
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cross sections of the radiative process e^e+ e^Peud'^ are calculated based 
on the exact amplitudes generated by GRACE and integrated numerically in 
five-body phase space using BASES. Again only the t-channel diagrams(non- 
annihilation diagrams) are taken into account. To compare the distributions 
the soft-photon corrections for the radiative process must be included. For 
this purpose QEDPS is implemented into the e'VeUdj calculations with a 
careful treatment to avoid a double-counting of the radiation effect. The 
definition of the hard photon is 

1. E^>1 GeV, 

2. opening angle between the photon and the nearest final-state charged 
particles is greater than 5°. 

The distributions of the hard photons are in good agreement as seen in Fig.|[ 
The total cross section of the hard photon emission agrees each other in 
2% level. We also calculated the radiative cross-section without soft-photon 
correction. If the soft-photon correction is not included in the radiative 
process, we find a 30% overestimation of the radiative cross-section with 
above experimental cuts. 

3.2 Cross sections with experimental cut 

It is also investigated how large the effects of the QED radiative corrections 
is for the single-VF production when realistic experimental conditions are 
imposed. The experimental cuts applied here are 6'e- < 5°, Mgg > 45 GeV, 
Efj_ > 20 GeV. For this study all the diagrams, not only the t-channel 
diagrams but also the s-channel, are taken into account in the calculations. 
The dominance of the signal from t-channel diagrams is 97%(90%) for the 
hadronic(leptonic) decay of the W boson, respectively. Total cross sections 
as a function of the CM energies at LEP2 with these cuts are shown in Fig.^ 
The QED radiative corrections on the total cross sections are found to be 7 
to 8% in this LEP2 energy range. If one uses the inappropriate energy scale 
s in SF, the ISR effect is overestimated by around 5%, which is larger than 
those of the no-cut case. SF with the proper energy scale shows a deviation 
from QEDPS around 0.5% for e'Veud process and 1.0% for e~z>e/Lt^f^. The 
agreement between QEDPS and SF becomes worse than the no-cut case, 
because the finite transverse-momentum of the emitted photons by QEDPS 
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changes the acceptance of the above electron- veto requirement. Hence for the 
reahstic experimental conditions the finite transverse-momentum of emitted 
photons should not be ignored. 

Finally a possible effect of the systematic error of the ISR effect on the 
anomalous TGC measurements is investigated. If the inappropriate energy 
scale is used in the ISR tools such as the CM-energy squared, the total cross 
sections with the experimental cuts includes 5% systematic error. Even if 
one used one of proper energy scales, there is 0.6% uncertainty on the cross 
sections due to ambiguity of the energy scale selection. These uncertainty of 
the total cross sections limits the experimental sensitivity of the anomalous 
TGC measurements. Total cross-sections of eveud and eVf-^iVy processes with 
experimental cuts as a function of the anomalous TGC is summarized in 
Fig.|^. The cross sections are obtained at the CM-energy of 200 GeV with all 
diagrams. A bound shown in dashed lines shows ±5% cross section variation 
from its standard model value. If one used the inappropriate energy scale, 
it affects ±0.05 on measurement and ±0.4 on A measurement. The 
ambiguity of the energy-scale selection gives the systematic error of less than 
0.01 on Aav and 0.1 on A. 



4 Conclusions 

The method to apply the QED radiative correction to the non-annihilation 
process was established. The conventional method, SF with the energy scale 
s gave about 4% overestimation in the LEP2 energies. The uncertainty due 
to the energy-scale determination was estimated to be about 0.6%. This 
uncertainty may affect the anomalous TGC measurements less than 0.01 on 
Ak and 0.1 on A. If one wants to look at the hard photon spectrum, the 
soft-photon correction to these radiative processes are needed. 

In this report we have treated two processes e~ e~Ueud and e~e^ 
e~z>e/i+z-'^. The CP conjugate processes e~e+ e'^v^ud and e~e~^ — > e+z/g/U"//^ 
give the same results and the other channels e~e+ — evecs and e~e^ 
eV(.TVr will show slightly different results due to their masses. On the other 
hand the self CP-conjugate process e~e^ —>■ e~z>ee^fe has an additional com- 
plexity, because two energy scales t_ = {pe- —qe-Y ^"^^ ^- = {Ve- —Qe- —Qd^Y 
can appear simultaneously. 
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Figure 1: Feynman diagrams of the process e e+ ^ e Ueud. First ten 
diagrams show the s-channel diagrams and last ten the t-channel ones. 
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Figure 2: Differential cross section of tlie total energy of emitted photon(s) 
obtained from the QEDPS (histogram) and from the SF(circle). 
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Figure 3: Total cross sections and those normalized by the QEDPS results 
for eveud and eVefJ^i^u processes without experimental cuts. The SF(t) denotes 
the SF with proper energy scale and the SF(s) with the inappropriate energy 
scale (s). Only t-channel diagrams(non-annihilation diagrams) are taken into 
account. 
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Figure 4: Differential cross sections of tfie liard plioton; Energy, transverse 
energy w.r.t. tfie beam axis, cosine of tlie polar angle, and opening angle be- 
tween photon and nearest charged-fermion. A histogram shows the QEDPS 
result and stars from the matrix element with soft-photon correction. Only 
i-channel diagrams(non-annihilation diagrams) are taken into account. 
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Figure 5: Total cross sections and those normalized by the QEDPS results 
for eueud and eue^iy^ processes with experimental cuts. The SF(t) denotes 
the SF with proper energy scale and the SF(s) with the inappropriate energy 
scale (s). All diagrams are taken into account. 
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Figure 6: Total cross-sections of eUeud and eUefii^u processes with exper- 
imental cuts as a function of the anomalous TGC. Dashed (dotted) Unes 
show ±5% (±0.6%) cross section variation from its standard model value. 
The cross sections are obtained at the CM-energy of 200 GeV with all dia- 
grams. 
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